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We present results for Bk calculated using HYP-smeared improved staggered fermions on the 
MILC asqtad lattices. In this report, the data is analyzed using the results of SU(2) staggered 
chiral perturbation theory (SChPT). We outline the derivation of the NLO SU(2) SChPT result, 
explain our fitting procedure, and outline how we estimate systematic errors. We also show the 
light sea-quark mass and lattice spacing dependence for both SU(2) and SU(3)-based analyses. 
Our preliminary result from the SU(2) analysis is B^NDR,^ =2 GeV) =0.512 ±0.014 ± 0.034 
and Bk = fiif(RGI) = 0.701 ±0.019 ±0.047. This is somewhat more accurate than our result from 
the SU(3) analysis. It is consistent with results obtained using valence domain-wall fermions. 
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1. Introduction 

This paper is the second in a series of the four proceedings describing our calculation of 
Bk- In the first, we provided a brief phenomenological introduction, outlined the results of SU(3) 
staggered chiral perturbation theory (SChPT), and explained the corresponding fitting strategy [j|]. 
Here, we explain briefly how we obtain the SU(2) SChPT result, outline how we use this result 
to analyze our data, and present a preliminary value, and error budget, for Bk- The details of the 
lattice ensembles and quark masses are as in Ref. [JXj] . 

The use of SU(2) ChPT was pioneered in the lattice context by the RBC collaboration [0]. 
One treats kaons and r\ 's as heavy, static sources for pseudo-Goldstone bosons (PGBs) composed 
of light quarks (the "pions"). Unlike the SU(3) version, SU(2) ChPT does not require an expansion 
in the ratio r s = wi^/Aqcd- The expansion parameter is thus r\ = ra^/AQCD (with mg the common up 
and down quark mass in our isospin-symmetric simulations). This improves convergence properties 
(as long as mi is light enough), although this comes at a price: one fits to a smaller number of 
data-points, and must do so with low-energy coefficients (LECs) that have an unknown, analytic 
dependence on m s j Kqcd. 

Our calculations have valence quark masses ranging from approximately m!? hys /10 to m\ ys . 
Despite the large upper value, we find that SU(3) fits give a reasonable description of our data. 
The main problem is the presence of multiple contributions to the NLO formula arising from dis- 
cretization errors or from truncation of the perturbative matching factors. These lattice artefacts 
are poorly determined in the fits and lead to ~ 5% uncertainties upon extrapolation to the physical 
light-quark mass. A major advantage of the SU(2) approach in our context is that all these artefacts 
are moved to NNLO, i.e. are known to be very small. This allows a better-controlled extrapolation 
to the physical mass. 

2. SU(2) Staggered ChPT Analysis 

The following description is necessarily very brief. A full explanation will appear in Ref. [||]. 

We need to generalize the SU(2) result for Bk given in continuum ChPT in Ref. [Q] to include 
the artefacts associated with staggered fermions. This generalization has been done in the SU(3) 
case in Ref. [Q], but requires a rather involved operator enumeration. 1 Rather than carry out a 
direct generalization to SU(2) SChPT, we instead have worked to all orders in r s (though to NLO 
in r£) within SU(3) SChPT. Using power-counting arguments, we find a simple result [J3j] : one can 
obtain the NLO SU(2) SChPT result simply by taking the SU(2) limit of the NLO SU(3) SChPT 
result, and then allowing (almost) all of the LECs to have an unknown dependence on r s . The only 
exception to this arbitrary dependence is the factor of multiplying the pion chiral logarithms, 
which remains unchanged. 

Applying this result we find the simplification noted above, namely that all the non-analytic 
contributions multiplied by LECs proportional to a 1 or a 2 are pushed to NNLO in SU(2) SChPT. 
This happens because the factor of M\ in the numerator of Bk is balanced by a chiral logarithm 
proportional to M\\o%(M n ). In SU(3) ChPT this ratio is of 0(1), but in the SU(2) case it is of 
NLO. The overall factor of a or a then moves the term to NNLO. 

'Our use of a mixed action leads to small corrections to the analysis of Ref. [0] that we have determined o,lllB||. 
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The final result is that taste-breaking effects only enter the NLO expression through the masses 
of the pions of different tastes. The predicted form is 

f th = d Q F Q + d\ + d 2 ^j + dy^- (2. 1) 

A Z X X 

where the chiral logarithms (defined in Ref. [jl]]) reside in the function 

Fo = 1 + 3^272{^ I ) + (L I -X I K~(X I )-2(£(X B ))} (2.2) 



<W> = ^ 2 \£(X l )+£(X P )+4£{X Y )+4£(X A ) + 6£(X T ) 



16 L- - - J ' (2 " 3) 

Here Xq (Lb) is the squared mass of the valence (sea) pion with taste B, which we know from our 
simulations or those of the MILC collaboration. The scale A x is arbitrary and we take it to be 1 
GeV. The coefficients dj have an unknown dependence on r s , and, in addition, at NLO do depends 
also on a 2 and a 2 . 

As for SU(3) fitting, we include a single analytic NNLO term — that proportional to d 2 - In the 
SU(2) case, however, we find that we can drop this term if we consider only the smallest valence 
light-quark masses. In the following we label such fits as "NLO", while if we include the d 2 term 
we label the fits with "NNLO". 

3. SU(2) SChPT Fitting 

A major advantage of the SU(2) analysis is that the fitting function is much simpler than that 
for SU(3). As a consequence, we do not need to include Bayesian priors for any of the parameters. 

Our choice of which valence quark masses to include in the analysis is exemplified by our 
approach on the coarse MILC lattices (a « 0.12 fm). Here the physical strange quark mass is 
am phy = 0.052, and we choose the nearest three values for the valence strange-quark mass, am y = 
{0.05,0.045,0.04}, in order to extrapolate to the physical value. For the valence down-quark mass 
we use our lightest four values, am x = {0.005,0.01,0.015,0.02}, to extrapolate to am p d ys . These 
choices ensure that the expansion parameter of SU(2) ChPT is relatively small: m x /m y < 1/2. 
Analogous choices of quark masses are made on the fine and superfine ensembles [||]. We call this 
choice "4X3Y", and use it for our central value. We have also considered "5X3Y" fits. 

In Fig. [lj, we show examples of the resulting fits. In the left plot, the (blue) octagons show the 



one-loop matched lattice data, which are then fit to the NNLO form of eq. (2.1). There are only 3 



parameters because Lp is fixed in this partially-quenched fit. With the fit parameters determined, 



we then evaluate the expression (gjj) at the physical pion mass and with taste splittings set to zero, 
i.e with Xq = L\ = M 2 . This removes all taste-breaking discretization and truncation errors, and 
results in the point shown with the (red) diamond. We call this procedure the "X-fit", and we repeat 
it for each of the values of am y . 

We then proceed to the "Y-fit", illustrated in the right panel of Fig. [lj Here we use either a 
linear or a quadratic fit to extrapolate the short distance to the physical value Fp = 2Af| —M 2 . We 
find, as seen in the figure, that the dependence is weak and close to linear. Thus we use the linear 
fit for our central value and the quadratic fit to estimate a systematic error. In the figure, the red 
diamond shows the result of linear extrapolation. 
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Figure 1: One-loop matched Bk fitted versus Xp (left) and Yp (right), on the MILC coarse lattices with 
amg = 0.01 and am s = 0.05. In the left panel am y = 0.05, and the (red) diamond shows the result after 
removing lattice artefacts. The fit type is 4X3Y-NNLO. Bk is obtained using one-loop matching. See text 
for more details. 



4. Dependence on the Light Sea Quark Mass 

After the X- and Y-fits, the resulting intermediate value for Bk still has the analytic depen- 
dence on the light sea-quark mass ami (entering through the term), as well as taste-conserving 
discretization and truncation errors. Here, we investigate the dependence on ami oc Lp, which we 
have studied on the MILC coarse ensembles. We also discuss the corresponding dependence of the 
result of the SU(3) SChPT fit discussed in Ref. 

In Fig. ||, we show the am^ dependence for both SU(3) and SU(2) fits. We have fit to the 
expected linear dependence [resulting in the (red) squares], and also, since the data show very little 
dependence on am^, to a constant [(red) diamonds]. 

The weakness of the dependence on ami 1S striking. It has the important consequence that our 
use of only a single light sea-quark mass on the fine and superfine lattices does not introduce a 
large uncertainty. We use the difference between the constant and linear fits as an estimate of the 
systematic error arising from the uncertainty in the ami dependence. As the figure shows, this error 
is somewhat smaller for SU(2) fitting than for SU(3). 

5. Continuum Extrapolation 

The dominant errors remaining at this stage are those due to taste-conserving discretization 
and truncation errors. These vary as a 2 a", where n = 0, 1 , . . . (n = is allowed since we do not use 
Symanzik-improved operators), and as a 2 . We cannot disentangle these effects using a fit to three 
lattice spacings, so proceed as follows. We fit our data to a linear function of a , and estimate the 
0{a 2 ) truncation error separately (as described in a companion proceeding [J7|]). More precisely, 
we fit to the data from the lattices with a = 0.12,0.09,0.06 fm having amj/am s fixed to 1/5. The 
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Figure 2: Bk (after X- and Y-fitting) versus ami for the SU(3) analysis (left) and for the SU(2) analysis 
(right). Data is from the MILC coarse lattices. The fit types are N-BT7 for SU(3) and 4X3Y-NNLO for 
SU(2). See text for details. 
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Figure 3: One-loop matched Bk plotted versus a 2 for the SU(3) analysis (left) and for the SU(2) analysis 
(right). The fit types are N-BT7 for SU(3) and 4X3Y-NNLO for SU(2). 



results for both SU(3) and SU(2) analyses are shown in Fig. [3| We use the extrapolated values for 
our main result, and take the difference between them and the results on the superfine lattices as an 
estimate of the systematic error due to the continuum extrapolation. 

The dependence on a 2 is noticeable for both analyses, with Bk increasing by about 6% and 
10% between a = and 0.12 fm in the SU(3) and SU(2) cases, respectively. Assuming a form 
gcpntj^ _|_ (aA) 2 ] this corresponds to scales of A as 400 and 500 MeV. These are reasonable values, 
indicating that HYP-smearing has reduced the discretization errors with staggered fermions down 
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to canonical size. 

6. Error Budget and Conclusion 



cause 


error (%) 


description 


statistics 


2.8 


4X3Y-NNLO fit (ensemble C3) 


discretization 


1.4 


diff. of (SI) and a = 


fitting (1) 


0.15 


X-fit: NLO vs. NNLO 


fitting (2) 


0.5 


Y-fit: linear vs. quadratic 


fitting (3) 


0.25 


constant vs linear ami dependence 


finite volume 


0.89 


20 3 (C3) versus 28 3 (C3-2) 


matching factor 


6.4 


AB^' (SI) 


n 


0.09 


uncertainty in r\ 



Table 1: Preliminary error budget forB^ obtained using SU(2) SChPT fitting. 



In Table [j], we summarize our present best estimates of the uncertainties in Bk arising from 
various sources. The method by which we estimate these errors is outlined in the "description", 
and has in most cases been explained above. The statistical error is obtained from using a global 
jackknife procedure. The finite volume estimate is obtained by comparing results on two volumes, 
as described in one of the companion proceedings [§]. This error is comparable to that estimated 
using NLO ChPT The error due to the use of one-loop matching is estimated in another of the 
companion proceedings The error due to the uncertainty in the scale r\ is estimated by varying 
the input values within the quoted errors and repeating the analysis. 

One error not included in this budget is that due to the strange sea-quark mass differing slightly 
from its physical value. Given the weak dependence on the light sea-quark mass, we expect this to 
be a negligible effect, but plan to make a more quantitative estimate in the future. 

Combining these errors, our current, preliminary result for Bk using SU(2) SChPT fitting is 

fi^NDR,^ = 2 GeV) = 0.512 ±0.014 ±0.034 [S U (2), PRELIMINARY] , 

B K = B K (RGl) =0.701 ±0.019 ±0.047 [S U (2), PRELIMINARY] , 1 ' 

where the first error is statistical and the second systematic. The total error is thus about 7%, and 
is dominated by the error in the matching factor. The total error is smaller than the 9% preliminary 
error we find for the SU(3) analysis [jl]]. The statistical error in the latter analysis is smaller (as can 
be seen in the earlier plots, and arises because there are more data points in the fits) but this is more 
than balanced by an increase in the systematic errors. 

Our result is consistent with other 2 + 1 -flavor unquenched results obtained using differ- 
ent fermion discretizations, all of which are summarized in Ref. [Js[] . For example the calcula- 
tion of Ref. [^] using domain-wall valence fermions on the coarse and fine MILC lattices finds 
B/-r(NDR,/i = 2 GeV) =0.527 (6) (20). The agreement between calculations using different fermion 
discretizations provides a highly non-trivial cross-check on the results. 

During the next year, we aim to reduce our errors by adding a fourth ("ultrafine") lattice 
spacing (a m 0.045 fm), adding other sea-quark masses, and by increasing the statistics. We are 
also working on two-loop and non-perturbative calculations of the matching factors. 
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